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1 Abstract 



The modulation of high-energy transients' (or steadily emitting sources') light curves due 
to the imperfect alignment of the detector's view axis with the spin axis in a spin-stabilized 
satellite is derived. It is shown how the orientation of the detector's view axis with respect 
to the satellite's spin axis may be estimated using observed light curves. The effects of 
statistical fluctuations are considered. 

Conversely, it is shown how the attitude of a spin-axis stabilized satellite as well as 
' the unknown position of a celestial source of high-energy photons may be determined 

using a detector whose view-axis is intentionally kept inclined and is known accurately 
beforehand. The case of three-axes stabilized satellites is also discussed. 
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2 Introduction: 

X 

In high-energy (X-ray and gamma ray) astronomy experiments photon detectors (scintil- 
lators such as Sodium Iodide, Cesium Iodide or gas-filled multi-wire proportional counters 
(MWPCs) or even solid-state detectors such as Si(PIN) or CdTe detectors) are flown on- 
board satellites. The satellites are either spin-axis stabilized (such as the Indian SROSS 
C-2 satellite) or three-axes stabilized. In the case of spin-axis stabilized satellites, nor- 
mally, the detector's view axis is aligned with the spin-axis of the satellite. However, 
sometimes this alignment is not perfect and the view axis of the detector makes a small 
but finite angle with the spin-axis of the satellite. In this case the observed light curve of 
a given celestial source of photons is modulated by the spin -period of the satellite. The 
amplitude of the modulation is a function of the polar angle (9 V ) between the detector's 
view axis and the satellite's spin axis. It also depends on the relative direction of the 
source with respect to the spin axis (9 ). 
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The spin-period is usually not strictly constant over the entire life time of the satellite. 
But over a small time period (during the entire duration of a high energy transient) it 
may be assumed to be almost constant. 

In the present paper we derive the amplitude of the light curve's modulation as a 
function of the space angle 9 V . We take a high energy transient light curve and superimpose 
the modulation (using a given value of V ). We show how the (if unknown) orientation of 
the view axis may be estimated reasonably accurately using the observed modulated light 
curve and the unmodulated light curve of the same high energy transient (say, a Solar 
X-ray flare) detected by another satellite. 

Finally, we discuss how this apparently disadvantageous and undesirable phenomenon 
may be useful in determining the position of a celestial high energy photon source given 
that the orientation of the detector's view axis (intentionally kept inclined) is known be- 
forehand. Also it is described how this same phenomenon might be utilised in determining 
the (unknown) attitude of the spacecraft by placing two small detectors (one inclined to 
the spin-axis, the other aligned with the spin-axis). 

3 Mathematical Formulation: 

Let XYZ (Fig.l) denote the celestial equatorial co-ordinate system. Let a s and 5 S denote 
the orientation of the satellite's spin axis in the XYZ system. We define 9 S = n — 5s and 
05 = cxs- Let the source direction be «o and Sq respectively. 




Figure 1: XYZ is the celestial polar co-ordinate system. OS is the spin-axis of the satellite. 
00' is the vector in the source direction. OV is the detector's view axis. The detector's 
view axis and the source rotate around the spin axis at the spin period of the satellite. 

Then, 9 = n — 5o and 0o = «o- The direction cosines of the source with respect to 
the spin-axis system may be obtained as 

sin9' cos(f>' \ I cos6scos(j)s cosOssin^s —sin9s \ ( sinOoCoscpo \ 
sin9' o sin(f)' = —sinfis cosfis sin9 siri(f) (1) 

cos9' J \ sin9scos(j)s sinOssinc^s cosQs J \ cos9 J 

(see, for example, Arfken and Weber). Hence the unit vector in the source direction is 

S = i sin9 cos(p + j sin9 siri(p + k cos9 (2) 
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Let the direction of the detector's view axis be 9 V and <p v respectively with respect to the 
spin-axis co-ordinate system. Then the unit vector along the view axis in the spin-axis 
system is given by 

V — i sin9 v cos(p v + j sin9 v sin<p v + k cos9 v (3) 
If we assume <p v — 0, this becomes 

— * — if —*/ 

V = i sin9 v + k cos9 v (4) 

This essentially implies that the spin-axis system is given a trivial rotation by an amount 
<f> v in the opposite sense. 

3.1 Determination of the inclination of the view axis: 

The modulation factor (assuming <f> v = 0) is given by 

M = sin9 v sin9 cos(ut + ip) + cos9 v cos9 (5) 

where we have put <p' = (uit+ip), ip being the epoch, t the time and u the angular velocity 
due to the spin. 

The modulation affects only the signal and not the background counts. The modulated 
counts 

Di = CiM = Ci(sin9 v sin9 Q cos(uti + ip) + cos9 v cos9 ) (6) 

where the Cj's are the counts in the unmodulated light-curve. 

To determine 9 V (9 is given) one has to consider another unmodulated light curve 
(this light curve has to be multiplied by an appropriate factor g such that the maxima 
(peaks) are the same for the two light curves. 

Therefore, we have 

sin9 v sin9 cos(ujti + ip) + cos9 v cos9 = fi (7) 

where fi = Di/Ci,Ci being known from a detector which is aligned (on-board a different 
satellite). Here, of course, the modulated light curve has to be multiplied by a suitable 
factor such that the maxima of the two light curves are equal. This may be written as 
(since 9 Q is known, i.e. both sin9 and cos9 are known). 

aixxcosiuti + + hx 2 = fi (8) 

where x\ = sin9 V) X2 = cos9 v . Also, a\ = sin9' , bi = cos9 . Taking the time average of 
both sides (over an integral number of cycles), 

< aiXicos(uti + ip) > + < bix 2 >=< fi > (9) 

aiXi < cos(uti + ip) > +bxx 2 =< f > (10) 

Since < cos(uti + tp) > is equal to zero for an integral number of cycles, the first term 
becomes equal to zero. Hence, 

hx 2 =< fi > (11) 

or, 

X2 = (l/6i) < fi > (12) 
Since, x 2 = cos9 v , 9 V = cos~ 1 (l/6 i ) < fi >. 
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3.2 Determination of Source Location 

If 9 V is known beforehand, it is possible to determine the location of a source in the sky 
using two detectors, one aligned, the other inclined (at a known angle) with the spin axis 
of the satellite. 

The equation for the modulated light curve (counts vs. time) is 

Di = Ci(sin6 v cos<j) v sin8 cos<f) + sin9 v sin<j) v sin9 sin<f)Q + cos8 cos8 v ) (13) 

where Cj is the counts detected during the ith time bin (when the view axis is aligned 
with the spin axis, i.e. 9 V — 0). The other symbols have their usual meanings. Assuming 
4> v = (reorientation of the spin-axis system), this reduces to 

Di = Ci(sin6 v sin6 cos<f) Q + cos9 cos9 v ) (14) 

Here, 4> Q = cos(uoti + ■0), ip being the epoch. For a detector whose view axis is aligned 
with the spin axis, 9 V — 0. Then 

A = C t (cos6' ) (15) 
From the last equation, for a given (say, the ith) time bin, 

cos9' = Di/d (16) 

which gives 

e' 0i = cos- 1 (D i /(C i )) (17) 



Let a = cos0' o (known) and sin9' = ±y (1 — a 2 )- If, for a second detector (for which 
9 V 7^ 0, i.e. sin9 v ^ 0, then 

D'^CKbr^ + ac) (18) 

This gives 

4 = co S - 1 ((l/&r)(/j;/C<;-ac)) (19) 

Here b = sin9 V) r = sin9' v and c = cos9 v . Knowing the values of 9' and 0q, the inverse 
transformation (corresponding to that given in eqn. (1)) may be used in order to obtain 
the values of 9 and O - 



/ sin9ocos(f>o 
sin9osin(f>o 
\ cos9q 



I cos9scos(j)s —siri(ps sin9scos(f>s 
cos9ssin(f>s cos(f>s sinOssin^s 
\ —sin9s cos9s 



I sin9' cos(j)' 

sin9' sin(j)' | (20) 
\ cos9' 



As evident there will be many such values (one set for each time bin) and their means 
and errors may be determined. 



3.3 Determination of Spacecraft Attitude 

If 9 V is known, it is possible to determine the attitude of the spacecraft as follows. This 
also requires two detectors-one aligned and the other inclined (at a known angle) with 
the satellite's spin axis. 
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The equation connecting the unknown direction of the spin axis, (9s and (fis) to the 
known (true) direction cosines and known (apparent) direction cosines of the source is 
the following. 

^ sin9' cos(f)' \ / cosOscoscps cos9ssin4>s —sin9s \ ( sin9 cos<f)o \ 

sin9' sin(f)' = —sirups cos4>s j sin9osiruf>o (21) 

\ cos9' J \ sin9scos(j)s sin9ssin4>s cos9s J \ cos9 J 

Three equations are obtained from the above matrix equation. One of the equations is 

— lisinips + hcoscps = k 2 (22) 

where l\ = sin9 cos(f) , l 2 = sin9 siruj) , and k 2 = sin9' siruj)' . This leads to the quadratic 
equation 

(l\ + l\)x 2 + 2k 2 l x x + (kl -ll)=0 (23) 

Here x = sirups- Solution of the last equation gives two values of x from which the value 
of (f>s may be obtained as 4>Si — s^ _1 (pi) and <ps 2 — sm _1 (p 2 ) where pi and p 2 are the 
two roots of equation (21). 

The other two equations are, respectively 

licos9scos(ps + l2Cos9ssiri(fis — hsinOs = k\ (24) 

and 

l\sin9scos(f)s + l 2 sin9s sirups + hcos9s = k% (25) 

Here k\ = sin9' cos(f)' , l 3 = cos9 and k 3 = cos9' . The last two equations give rise to the 
following equation 

{k\ + kl)y 2 + 2hl 3 y + {l\ - kl) = (26) 

Here y = sin9s- Solutions of the above equation gives two values of y from which the 
values of 9s may be obtained as 9s 1 = sm _1 (gi) and 9s 2 = sm~ 1 (g 2 ) where qi = and q 2 = 
are the solutions of eqn. (24). 

The correct values of 9 S and <ps are to be decided based on physical considerations. 

3.4 Three Axes Stabilized Spacecrafts 

In the case of a three axes stabilized spacecraft one may have a small spinning platform on 
an extended boom on which the aligned and inclined detectors may be placed. The above 
mentioned procedures may be utilised to determine either source location or spacecraft 
attitude in this case as well. 

4 Results: 

4.1 The Timing Data and the Light Curve: 

Due to the non-availability of real data at hand, we take recourse to generating data 
artificially. 

The total duration of the time series data (detected counts vs. time) is equal to 512 
seconds with a time resolution of 256 ms (this is the width of each time bin). 
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During the leading 51.2 seconds and the trailing 51.2 seconds the detected counts are 
due only to the background. The background has a mean value of 1.8 counts and is 
fluctuated according to a Poisson distribution. 

Between the leading and trailing background data, in the remaining time interval of 
409.6 seconds a light curve is generated artificially. The light curve is essentially triangular 
in shape. The total duration (409.6 seconds) of the light curve is divided into two parts. 
The rising portion of the light curve has a duration of 153.6 seconds. The decaying 
portion has a duration of 256 seconds. Both the rising portion and the decaying portion 
of the light curve are linear in shape. The rising portion has a slope of +3.8 while the 
decaying portion has a slope of -2.28. Poissonian fluctuations are superimposed on the 
light curve generated by the Monte Carlo method (both in the rising as well as in the 
decaying portion). This simulated light curve is shown in Fig.2. 
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Figure 2: The unmodulated light curve. The leading 51.2 seconds and the trailing 51.2 
seconds consist of only the background counts with an average value of 1.8. Poissonian 
fluctuations are superimposed on the background counts. The intervening 409.6 seconds 
constitute the simulated light curve due to the high energy transient source. The rising 
portion has a duration of 153.6 seconds and has a slope Of 3.8. The decaying part has a 
duration of 256 seconds and has a slope of -2.28. 



4.2 The Modulated Light Curve: 

Due to the spinning motion of the satellite and the fact that the detector's view axis is not 
aligned with the satellite's spin axis, the light curve of the high energy transient (it is true 
also in the case of a steadily emitting celestial object) will be modulated in amplitude. 
The modulation factor is calculated as described earlier. 

The modulated light curve corresponding to the original (unmodulated) light curve 
(Fig.2) is shown in Fig. 3. Here 9 S is the polar angle and <ps is the azimuth angle of the 
source as seen from the spin-axis frame of reference. 6 V and <f> v are the corresponding 
parameters for the view axis. To simplify the problem we assume <p v = 0. This essentially 
means a trivial rotation of the spin axis coordinate system about its Z-axis by an amount 
equal to <f> v . 

The results of determination of 9 V are displayed in Table 1. The error on the mean of 
the estimated values of 9 V = 0.283 degrees which is equal to 1.25 percent. 
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MODULATED LIGHT CURVE 
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Figure 3: The simulated light curve in Fig.l of the high energy transient modulated by 
the spin period of the satellite. The spin period is equal to 11.78 seconds. The counts in 
the light curve are fluctuated according to a Poisson distribution. 



Table 1 



Table 1: Estimated values of 9 V and their errors. 
Range of Data True Value of Estimated value Error on 



(in cycles) 9 V of 9 V 9 V 

(in degrees) (in degrees) (in degrees) 

I- 2 22.6 22.607 0.007 

3-4 22.6 23.580 0.980 

5-6 22.6 23.479 0.879 

7-8 22.6 22.407 -0.193 

9-10 22.6 21.600 -1.000 

II- 12 22.6 22.231 -0.369 
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5 Discussions: 



Obviously the amplitude of the modulation of a high energy transient (or even a steady 
source) light curve depends on the parameters 9 V , 9$ and <p v etc (as shown in the figures 2 
and 3). Since one assumes that the average of the term cosut, < cosut > equals zero, one 
has to consider the time-history (light-curve) data only for an integral number of the spin 
period of the satellite. In the present work 2 cycles of data has been used for estimating 
each value of 6 V . 

It has been shown that even in the presence of fluctuations, this method is still viable 
and yields quite good results. 

However, the accuracy in the estimated value of 9 V depends on the background. In 
another calculation where a background has been assumed which is 13 times larger, the 
errors in the estimated value of 9 V also becomes much larger, typically a few degrees. 
Therefore, this method may be used effectively in only those cases where the signal-to- 
noise ratios are quite large. 

The spin period of the satellite and the integration time of the light curve are very 
much realistic. Actually, the GRBM (Gamma Ray Burst Monitor) on-board the SROSS 
C-2 satellite possessed an integration time of 256 ms (Sinha, S et al) although the spin- 
period was nearly 12 seconds (about half of the value used in the present work). 

In the present work one assumes small (ideally point size) detectors since the math- 
ematical analysis does not assume any finite value of A<fi (due to the finite size of the 
detector). It is quite fortunate that presently solid state (Si(PIN) and CdTe) detectors 
are available having size as small as lmmXlmm that may be used to detect X-ray flares 
from the Sun. 

6 Conclusions: 

The effects of modulation of light curves of celestial high energy photon sources when the 
detector's view axis is not aligned with the spin axis of the satellite has been described. 
It is shown how the orientation of the view axis with respect to the spin axis may be 
determined by comparing this modulated light curve with the light curve obtained from a 
detector whose view axis is perfectly aligned with the spin axis (either in the same satellite 
or in a different satellite). Finally, the usefulness of this apparently disadvantageous 
situation is described- how this effect may be utilised in order to estimate the attitude of 
the spin axis or to determine the unknown position coordinates of a celestial high energy 
photon source. Further work in this direction will be reported shortly. 
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